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ABSTRACT 

We report new observations of the 1 io — 1 i i (6 cm) and 2i\— 2yi (2 cm) transitions of ortho-formaldehyde 
(0-H2CO) in absorption at z = 0.68466 toward the gravitational lens system B02 18+357. Radiative transfer 
modeling indicates that both transitions are anti-inverted relative to the 4.6 K cosmic microwave background 
(CMB), regardless of the source covering factor, with excitation temperatures of ~1 K and 1.5-2 K for the 
1 10-1 1 1 and 2\\—2\2 lines, respectively. Using these observations and a large velocity gradient radiative transfer 
model that assumes a gradient of 1 kms" 1 pc -1 , we obtain a molecular hydrogen number density of 2 x 
10 3 cm" 3 < n(H 2 ) < 1 x 10 4 cm" 3 and a column density of 2.5 x 10 13 cm" 2 < JV(o-H 2 CO) < 8.9 x 10 13 cm" 2 , 
where the allowed ranges conservatively include the range of possible source covering factors in both lines. The 
measurements suggest that H2CO excitation in the absorbing clouds in the B02 18+357 lens is analogous to that 
in Galactic molecular clouds: it would show H2CO absorption against the CMB if it were not illuminated by 
the background quasar or if it were viewed from another direction. 

Subject headings: galaxies: individual (B02 18+357) — galaxies: ISM — radiation mechanisms: non-thermal 
— radio lines: galaxies — quasars: absorption lines 



1. INTRODUCTION 

Formaldehyde (H2CO) centimeter transitions have been ob- 
served with excitation temperatures less than the tempera- 
ture of the cos mic microwave bac kground (T ex < Tcmb) m 
the Milky Way (Pal mer et al.ll 1969b . The "anti-inversion," or 
non-thermal absorption, is pumped by collisions, coupling 
H2CO excitation to the number density n(H2) of molecu- 
lar hydrogen even at kine t ic temperatures as low as ^10 K 
(lTownes& Cheung! 119691 : lEvans et al.1 119751 : iGarrison et al] 
119751) . Thus. HtCO provi des a density and temperature di- 
agnostic for molecul ar gas dMiihle. Seaquist. & Henkel2007l : 
IMangum et al.l 12008). and it often produces detectable lines 
before dense molecular clouds begin to form stars and warm 
sufficiently to emit lines from the usual molecular tracers such 
as CO and HCO + . 

The rotation states of H2CO are defined by the three quan- 
tum numbers Jk„k c (total angular momentum and the angu- 
lar momenta about the axes with the smallest and largest mo- 
ments of inertia, respectively). The lower energy transitions 
for K a < 2 are shown in Figure Q] The spin symmetry of 
the hydrogen nuclei defines the ortho (K a = odd) and para 
(K a = even) states. Since the spin symmetry is about the a 
axis, which is aligned with the electric dipole along the C- 
O double bond and thus unchanged by dipole interactions, 
the mean transition time from 0-H2CO to P-H2CO states (or 
vice versa) is longer than the mean life of the molecule, and 
the ortho/para ratio preserves information about the domi- 
nant formation channel: a ratio of unity indicates catalysis 
on the surface of cold (Tki n < 15 K) dus t grains, while a ra- 
tio of 3 indicates gas-p hase formation dKahane et al.l 119841 : 
iDickens & Irvind[l999l) . 

The most commonly detected 0-H2CO lines are the AJ = 0, 
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AK C = ±1 "^-doublet" li -ln 4.8 GHz (6 cm) and 2 n -2 12 
14.5 GHz (2 cm) centimeter transitions and the AJ = ±1, 
AK C = ±1 millimeter transitions connecting the 7=1 states 
to the J = 2 states at 140.8 and 150.5 GHz (2i 2 -ln and 
2n— lio, respectively). The coupling of the centimeter and 
millimeter transitions allows the /^-doublet transitions to be 
driven significantly out of local thermodynamic equilibrium 
(LTE) and makes H2CO an in situ probe of conditions in 
molecular clouds. The re lative excitation of ^ -doublets is 
a powerful densitometer (IMangum et al.l 120081) . while the 
relative e xcitation of millimeter lines is a molecular ther- 
mometer dMiihle. Seaquist. & Henkell 120071) . The 1 10 — 1 1 1 
and 2\\—2\2 transitions have both been observed with pop- 
ulation anti-inversions, causing stimulated absorption of the 
CMB, and the lip— In transition has been observed as a mase r 
(lHaar & Pellingl[l97l lEvans et aU [19751 : lArava et al.ll2008h . 

We report observations of the lio— In and 2\\—2\<x transi- 
tions in absorption toward B02 18+357, a gravitational lens 
system with two images sepa rated by 0"33. The lens lies at 
a redshift of z = 0.68466(4) dBrowne et alJI 19931) while the 
background quasar is at z = 0.944(2) (ICohen et al 1120031) . At 
low frequencies, the "gravitational lens with the smallest sep- 
aration" of its images also displays the smallest known Ein - 
stein ring, with a diameter of 0"335 dPatnaik et all 119931) . 
The two images, dubbed A and B in order of radio bright- 
ness, display significantly different absorption characteristics, 
with B unobscured while A shows strong molecular absorp- 
tion lin es and is heavily obscured by dust in optical obser - 
vations dWiklind & Combes|[l995l: iGrundahl & Hjorthll 19951) . 
Both A and B show the background source to have a core-jet 
morphology. 

Many molecules have been observed at the redshift of the 
lens in absorption against the strong continuum of image A, 
including CO, HCO + , HCN, H 2 0, NH 3 , OH, and H 2 CO; typ- 
ical line widths are 10-15 km s" 1 dWiklind & Combeslll99l 
Combes & Wlklrndl \l99% iHenkel et al J 120051 : iKanekar et aU 
20031 : iMenten & Reidl 1 19961) . H I observations indicate 
the lensing galaxy to be gas-rich, with a column density 
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Table 1 

Journal of o+UCO Observations Toward B0218+357 



Telescope 


Receiver 


UT Date 


Transition 


Vo 


fobs 


t ira 


Av a 


rms 










(GHz) 


(GHz) 


(s) 


(kms" 1 ) 


(mJy) 


Arecibo 


S-low 


2004 Aug 3-4 


1 10 — 1 1 1 


4.829660(1) 


2.8668455(6) 


3600 b 


1.3 


0.8 


GBT 


X 


2006 Apr 26 


2n-2i2 


14.488479(1) 


8.6002392(6) 


3280 


2.1 


1.1 



" Av is the rest-frame spectral velocity resolution of the final, smoothed spectrum. 
b Two 1 800 s sessions from consecutive days were combined for the final spectrum. 
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Figure 1. Formaldehyde (H2CO) energy level diagram showing only the 
lower levels and K a < 2. The indicated transitions are in GHz. 



10 18 (7;//) cm" 2 across a line width (FWHM) 
where T s is the spin tem perature and / is the 



N(U I) = 4 x 
of 43kms~\ 

H I covering factor (ICarilli et al.l 17993). Four spectrally re 
solved absorption components can be identified in front of im- 
age A in high-reso lution interferometric HCO + observations 
(lMulleretal.1120071) . 

H 2 C O was first observed in B02 18+357 bv lMenten & Reid l 
(1996) in the 2n-2i2 transition in absorption. iJethava et alJ 
(120071) followed these observations with a six-line study of 
H2CO in B02 18+357, adding observations of the lowest two 
A7 = 1 rotational transitions of both the ortho and the para 
species, as well as the 1 10 — 1 1 1 transition of 0-H2CO. Us- 
ing th e 55 K kinetic tem perature derived from NH3 observa- 
tions dHenk el et al. 2005) and a large velocity gradient (LVG) 
model to determine the non-thermal H2CO excitation due to 



collisions with H2, IJethava et al.l ( 120071) determined a spatial 
number density of n(H2) < 10 3 cm -3 and a column density 
of A^rUCO) « 5.2 x 10 13 cm" 2 (no confidence interval is re- 
ported) in the strongest absorption component, with a weaker 
component adding ~15% to the column density. They mea- 
sured an o-H 2 CO/p-H 2 CO ratio of 2.0-3.0 with a best fit of 
2.8, giving Af(o-H 2 CO) = 3.8 x 10 13 cm -2 . They do not ad- 
dress the issue of H2CO excitation temperatures, which are of 



interest because anti-inversion of the lines would imply that 
the gas in the B02 18+357 lens would be visible in absorption 
against the CMB even without a background quasar to illumi- 
nate it. Obtaining excitation temperatures of H2CO lines in 
B02 18+357 will help to set H2CO observability conditions in 
non-le nsing (or otherwise un illuminated) galaxies. 

The lMangum et all (120081) survey of local galaxies demon- 
strates the feasibility of determining both «(H2) and 
Af(o-H2CO) strictly from the lio— In an d ^iv^-u transitions. 
The relative strengths of the transitions are set by the degree 
of collisional excitation, giving a precise densitometer that, 
in turn, yields the excitation temperatures of the transitions. 
In our analysis, we use the two-line densitometry method 
with new data, including 1 io~l 11 observations of higher spec- 
tral resolution and s ignal to noise ratio than that obtained by 
IJethava et al.l (120071) and 2n-2!2 observations with lower sig- 
nal to n oise (4.7 versus 8.8) b ut comparable resolution to 
those of iMenten & Reidl (119961) . Our analysis, while also 
yielding A^o-ITCO) and n(H2) from an exhaustive search of 
the density space with our LVG model, focuses on line exci- 
tation temperatures in B02 18+357 to show that the o-JTCO 
observed in absorption is anti-inverted relative to the 4.6 K 
Tcmb at z = 0.68, making it analogous to Galactic molecular 
clouds that show anti-inversion of the /if-doublet transitions 
and thus CMB absorption. 

In Section [2] we describe our observational and data re- 
duction procedures used to acquire and reduce the new high 
resolution spectra of the 1 10 — 1 1 1 and 2n-2!2 0-H2CO transi- 
tions in B02 18+357. The properties of both lines are reported 
in Section [3] and we use those results to determine «(H-2), 
Af(o-H2CO), and line excitation temperatures in Section [4] 
We show that the absorbing clouds in B02 18+357 are anal- 
ogous to Galactic molecular clouds that show anti-inversion 
of 0-H2CO centimeter lines relative to the z = Tcmb> and 
that the anti-inversion of /if-doublet transitions persists up to 
the 615-616 (101 GHz rest) transition (Section|5]). The conse- 
quences of our results are discussed in Section [6] 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Arecibo Telescope: 1 10 — In Line 

We observed the lio-ln 4.829660(1) GHz transitior0of o- 
H 2 CO, redshifted to 2.8668455(6) GHz, toward B0218+357 
at the Arecibo radio telescope0in August 2004 (TableQ}- Ob- 
servations were conducted with a 6.25 MHz band centered on 
the redshifted line, using 5 minute position-switched scans 
with a calibration diode fired after each position-switched 
pair and spectral records recorded every 6 s. Total integra- 
tion time was 3600 s. The autocorrelation spectrometer used 

5 All line freq uencies are from the JPL molecular line database 
IPickettet all 19981) , 

6 The Arecibo Observatory is part of the National Astronomy and Iono- 
sphere Center, which is operated by Cornell University under a cooperative 
agreement with the National Science Foundation. 



Formaldehyde Anti-Inversion in B02 18+357 



3 



Table 2 

Measured and Calculated Formaldehyde Line Properties 



Transition 




Compt 


Vobs-Vsys 


Depth 


FWHM 


Tapp 


JVappdv 




(GHz) 




(km s _1 ) 


(mJy) 


(km s -1 ) 




(kms -1 ) 


1 10 _ 1 11 


4.829660(1) 


1 


+3.7(0.4) 


9.6(0.4) 


12.3(0.6) 










2 


-9.8(2.4) 


1.8(0.3) 


14.2(4.0) 










Total 


+3.4(0.4) 


10.8(0.8) 


12.6(0.6) 


0.017(1) 


0.239(8) 


2n-2i2 


14.488479(1) 


1 


+4.2(0.8) 


5.2(1.1) 


8.5(1.5) 


0.008(2) 


0.083(10) 



Note. — v b s - v sys is the rest-frame centroid velocity offset from the assumed systemic redshift, 
z = 0.68466, FWHM is the rest-frame line width, and r app is the maximum apparent optical depth in the 
line, assuming a covering factor of unity and a continuum in image A of 650 ± 1 6 and 663 ± 22 mJy in 
the 1 io — lii and 2\2 lines, respectively (Section^). Uncertainties are strictly statistical uncertain- 
ties from the spectra and do not reflect the systematic uncertainties associated with flux calibration or 
the radio continuum determinations. Reported quantities for the lio-lii "Total" and lines are 

measured directly from the data and are not dependent on profile fitting; the individual components of 
the 1 iQ-1 1 1 line are from Gaussian fits. 



nine-level sampling in two (subsequently averaged) polariza- 
tions. Bandpasses were divided into 1024 channels and Han- 
ning smoothed to 512 channels. Rest-frame velocity resolu- 
tion was 1.3 km s" 1 (see Table [T]). The observed band was 
interference-free in the vicinity of the observed line. 

Records were individually calibrated and bandpasses were 
flattened using the calibration diode and the corresponding 
off-source records. Records and polarizations were subse- 
quently averaged, and a polynomial baseline with variations 
much larger than the line was fit to and subtracted from the 
spectrum. Systematic flux calibration errors in these data 
are of order 10%. All data reduction was performed in the 
Astronomical Information Processing System++, AIPS++Q 
Arecibo spectra often show spectral standing waves due to 
resonances of the strong continuum flux densities within the 
telescope superstructure. The line width was significantly less 
than the size of the standing wave features, so the spectrum 
was not significantly affected. 

2.2. Green Bank Telescope: Line 

We observed the 2 U -2 12 14.488479(1) GHz transition of 
H 2 CO, redshifted to 8.6002392(6) GHz, toward B0218+357 
with the Green Bank Telescope^ (GBT) on 2006 April 26 (Ta- 
ble|T}. Observations were conducted with a 50 MHz bandpass 
centered on the redshifted line in 5 minute position-switched 
scans with spectral data recorded every 3 s and a winking 
calibration diode firing during every other record. The to- 
tal on-source integration time was 3280 s. The autocorre- 
lation spectrometer used nine-level sampling in two (subse- 
quently averaged) linear polarizations. Bandpasses were di- 
vided into 16384 channels, Hanning smoothed to 8192 in- 
dependent channels, and 10-channel Gaussian smoothed to a 
rest-frame velocity resolution of 2. 1 km s -1 . Records were in- 
dividually calibrated and bandpasses were flattened using the 
calibration diode and the corresponding off-source records. 
The individual spectra were interference-free and fairly flat, 
but the DC levels of the individual scans showed significant 
fluctuations of order 10%. Scans and polarizations were sub- 
sequently averaged, and a fifth-order polynomial baseline was 
fit to the full 50 MHz bandwidth and subtracted. Systematic 
flux calibration errors in these data are of order 10%. All data 

7 AIPS++ is freely available for use under the Gnu Public License. Further 
information may be obtained from http://aips2.nrao.edu 

8 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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Figure 2. Formaldehyde lio — In (6 cm) and 2n — 2i2 (2 cm) absorption to- 
ward the gravitational lens system B02 18+357. The rest-frame velocity scale 
assumes a heliocentric redshift of z = 0.684660, and the spectral resolutions 
are 1.3 and 2.1 km s -1 in the lower and upper lines, respectively. The dashed 
line shows a two-component G aussian fit to the 6 cm li ne profile. The 2 cm 
line spectrum, also detected by Menten & Reid 1 1996), is offset by 4 mJy, 
with the zero point indicated by the dotted line. 

reduction was performed in GBTIDL0 
3. RESULTS 

The Arecibo spectrum of the 1 01 - In transition requires 
two Gaussian components for a good fit to the line profile 
(Figure 0. The rest frame velocity resolution is 1.3 km s -1 . 
Table [2] lists the two fit components and the properties of the 
total line. The total line properties are not based on Gaussian 
fits; they are computed directly from the spectrum. All veloci- 
ties are in the rest frame with respect to a heliocentric redshift 

9 GBTIDL (http : //gbtld l ■ nrao ■ edu/t is the data reduction pack- 
age produced by NRAO and written in the IDL language for the reduction of 
GBT data. 
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of z = 0.68466. The integrated optical depth in this transition 
is 0.239(8) km s -1 over -20 km s" 1 < v < +17 km s" 1 assum- 
ing a continu um o f 6 50 ± 16 mJy and a unity covering factor 
(see Sections l4~T1andl4~2l). 

The GBT detection of the 2\\ — 2yi transition has a lower 
signal to noise ratio than the Arecibo lio— ln transition. De- 
spite this and the lower spectral resolution of the GBT spec- 
trum (2.1 versus 1.3 km s" 1 ), both lines display similar pro- 
files (Figure |2j, although the 2\\—2\2 detection is insufficient 
to justify a second profile c omponent as s een in the 1 io — 1 i i 
spectrum and in the data of iJethava et al.l (120071) . The prop- 
erties listed in Table |2] for this line have been measured di- 
rectly from the data. The integrated optical depth, spanning 
the same velocity range as the 4.8 GHz line, is 0.083(10) km 
s" 1 , assuming a continuu m of 663 ±22 mJy and a unity cov- 
ering factor (see Sections [4. 1134.21 ). 

The observed profiles are narro wer than the 15 km s" 1 ob - 
served in CO, HCO + , and HCN dCombes & WiklindlfT996h . 
but broader than the NH3 lines that range from 6.6(1.6) to 
9.3(1.8) km s" 1 dHenkel et al.ll2005l) . OH spec tra (1.67 GHz 
rest) clearly show a two-component profile dKanekar et alJ 
120031) . although at 20 and 27 km s" 1 the line components are 
broader than the ~10 km s" 1 H2CO lines. Our detection of 
the 2 1 1 —2 17, transition i s com parable to that of iMenten & Reidl 
d!996l) . IJethava et al.l (120071) report a much narrower 1 io~l 1 1 
line — with a FWHM of 5.5(5) km s" 1 — but our high sig- 
nal to noise ratio and resolution make our measurement more 
statistically robust. Both of our observations are coinci- 
dent in velocity space with previous molecular detections in 
B02 18+357. 

4. ANALYSIS 
4.1. Source Continuum 

In general, while the continuum determinations may be un- 
certain and lead to errors in the absolute optical depth and to- 
tal column density, the computation of quantities that depend 
on ratios of (low) optical depths, such as excitation tempera- 
tures and n(H2), will be significantly less affected, provided 
that systematic errors are consistent across frequency. 

The continuum emission at 2.9 GHz (the frequency of the 
redshifted 4.8 GHz I10-I11 transition) is nearly flat as a func- 
tion of frequency for images A and B dO'Dea et al.l [1991; 
iPatnaik et al.l 1 1993b but steeper for the Einstein ring and the 
extended emission halo. At 5 GHz, the ratio of the two com- 
pact flat spectrum components A:B is ~ 3.0 and the Ein- 
stein ring emission accounts for 10%-20% of the total emis- 
sion dO'Dea et al.ll 1 992t IPatnaik et al.lll993l) . The two sources 
show time variability with a 10.5 day period, but the variabil- 
ity is constrained to be < 8% at low frequencies (Mi ttal etaT] 
2006). Extended emission is detected at 1.465 and 1.63 GHz 
at the ~ 10% level , and additional stru cture is seen at 5 GHz 
dO'Dea et al.lll992l : lPatnaik1^[l99ll) . 

Comparison to other molecules co nfirms the association 
of H2CO absorpti on with image A dMenten & Reidl 119961: 
iMuller et alj|2007l) . The appropriate continuum level to use 
in apparent optical depth calculations for the 1 10 — 1 1 1 transi- 
tion is thus the total continuum at 2.9 GHz (8.6 GHz for the 
2n-2i2 line) less image B, the Einstein ring emission, and 
any extended source structure. Interpolating the flux densities 
for images A and B and the extended emission from interfer- 
ometric observations allows measurement of the fraction of 
the continuum from image A, a fraction crucial to disentan- 
gling the absorption in front of image A from the unabsorbed 



continuum of image B and other extended structure. 

The flux density of image A at 2.9 GHz can be interpo- 
lated from a linear fit (in the log) to the flux density ver- 
sus frequency of image A. We use th e very long baselin e in- 
terferometry (VLB I) observations of iMittal etall d2006l) be- 
cause these straddle 2.9 GHz and because these observa- 
tions, of all the available interferometric measurements, were 
taken closest in time to our spectra (2002 January versus 
2004 August for the Arecibo 1 10 — 1 1 1 observations and 2006 
April for the GBT 2\\-2\n observations). The continuum, 
interpolated from a fit to the image A data points spanning 
2.25-15.35 GHz (the 1.65 GHz datu m is dropped due to a 
steep downturn below 2.25 GHz; see IMittal et al.l 120061) . is 
650 ± 16 mJy, which neglects systematic and calibration er- 
rors typically of order 10%. The image A continuum is nearly 
flat, and interpolatio n error is le s s of a concern than over- 
all flux calibration. M ittal et al.l d2006|) degrade their reso- 
lution to 50 x 50 m as, comparab l e in r esolution to MER- 
LIN observations by IPatnaik et al.l (119931) to account for the 
loss of flux in VLBI observations, yielding similar flux den- 
sities to within 10%. Note that extrapolations from lower 
resolution (and much older) observ a tions p redict higher flux 
densities at 2.9 GHz: lO'Dea et al.l d!992l) predict 868 ± 64 
mJy, and IPatnaik et al.l d!993l) predict 791 ± 68 mJy. We 
adopt 5a(2.867 GHz) = 650 ± 16 mJy, bearing in mind that 
this estimate could potentially have large systematic errors. 
For the 14 GHz line redshifted to 8.6 GHz, we use the 
same VLBI data and continuum estimation method, adopt- 
in g S A (8.600 GHz) = 663 ± 22 mJy. The VLBI observations 
of IMittal et"ai1 d2006l) are degraded to 10 mas at 8.4 GHz and 
5 mas at 15 GHz. 

4.2. Source Covering Factor 

Calculation of the true optical depth of each line requires 
knowledge of the covering factor /„ describing the fraction of 
image A obscured by the absorbing cloud at the line frequency 
v. The covering factor corrects the observed line intensity to 
account for dilution of the line due to incomplete coverage 
of the continuum source by the absorbing cloud. Thus, the 
optical depth T\ m<t of the line takes the form 

where S\ mt and S v are the line and continuum fluxes, re- 
spectively. N4iU2m^te£_obs^rvarions at v > 100 GHz show 
/„ > 0.77 dWiklind & Combeil 19951) . consistent with a cover- 
ing factor of unity in front of image A and in f ront of image 
B; the flux ratio (A/B) is close to 4 at 100 GHz dMuller et alJ 
120071) . 

Due to th e increased solid angle of image A at lower 
frequencies dMittal et al.ll2006| ). the covering factors for the 
1 10-1 11 and 2n-2i2 transitions could be significantly less 
than unity even a fter excluding fl ux fro m image B and ex- 
tended emission. Kane kar et alJ (J2003) estimate the cover- 
ing factor to be 0.4 at 990 MHz (redshifted 18 cm OH), 
and IJethava et alJ d2007l) estimate the covering factor to be 
of the form f u = (i^/lOO) 1 for v in GHz, determining that 
0.0 < x < 0.5 by interpolating between the millimeter obser- 
vations of IWiklind & Combe si d 19951) and the 990 MHz OH 
observations of iKanekar et al.l (120031) (Note that, in arriving 
at this range, they refer to the OH lines as being at 850 MHz, 
the redshifted freq uency of the H I line.) For discussion, 
IJethava et al.l d2007l) adopt the most extreme frequency depen- 
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dence of x = 0.5, which gives covering factors of 0. 17 and 0.29 
for the 1 io — 1 i i and 1\\—2yi lines; the results they report are 
determined significantly by this choice of covering factors. 
Since this paper is concerned with excitation temperatures, 
which are determined via the ratio of the two lines, properly 
accounting for the frequency dependence of the covering fac- 
tor is vital to obtaining accurate results. 

VLBI observations of lMittal et al.l (120061) show an extended 
source at low frequencies that splits image A into two point- 
like images separated at 15.4 GHz by ^1.4 mas, correspond- 
ing to a physical separation in the lensing galaxy of 9.6 pc 
(for a flat universe with Ho = 73 km s _I Mpc -1 , D, m = 0.27, 
and fi v = 0.73). The solid angle il A follows a power law of in- 
dex of a « -2.7 (where CIa oc v a ) as frequency increases from 
1.65 to 15.65 G Hz. (By comparing centimeter and millime- 
ter observations. iHenkel et al.l(l2005h argue that -2 < a < — 1.) 
The relatively flat spectrum of the south-western source in im- 
age A suggests that it is the flat-spectrum core of the quasar 
and the location of continuum emission seen in millimeter ob- 
servations, where the frequency -dependent jet would not be 
observable dPatnaik et al.ll 1995b . Since the core is fully ob- 
scured (f v = 1) in millimeter observations, where the jet is not 
visible, the absorbing cloud lies along — but is not necessarily 
centered upon — the line of sight to it. 

At 8.6 GHz, where the jet and the core are not distinctly sep- 
arated, the ~15 mas 2 r egion containing the core emits ~60% 
of the flux of image A dPatnaik et al.l d 19951) measure 62(1)% 
with VLBI observations at 15 GHz). The total area of im- 
age A with flux density 5 times the rms backgr ound, found 
by extrapolating from the VLBI observations of Mi ttal et alJ 
(120061) . is ~32 mas 2 at this frequency. A spherical molecular 
cloud of 30 pc diameter subtends ~15 mas 2 (~47% of 0^ at 
8.6 GHz) at z = 0.68. Such a cloud gives a minimum covering 
factor of 0.3 if the core, which has /„ = 1 at millimeter fre- 
quencies, is obscured only by a limb of the cloud; due to the 
comparable line wid ths observed at both centimeter and mil- 
limeter frequencies ( Jet hava et al.l 12007). it is likely that the 
absorbing cloud is fortuitously aligned with image A and that 
the covering factor is close to 0.6 for a 30 pc cloud. This is a 
strictly geometric argument, and weighting the covering fac- 
tor by the fraction of flux that is obscured will give /145 > 0.6 
even for somewhat smaller clouds. 

The diameter of the cloud is likely to be significantly 
larger than 30 pc, since the optical obscuration of image 
A suggests the presence of large amounts of dust and pos - 
sibly a giant molecular cloud dGrundahl & Hjorthl 119951) ; 
IWiklind & Combesl d 1999T) report a visual extinction of A v = 
850 in front of image A. Applying the relation between cloud 
size S (equivalent to diameter for a spherical cloud) in pc and 
velocity dispersion a v in km s" 1 measured e mpirically with 
CO for Galactic clouds (ISolomon et al.ll 19871) of 

/ t\ \ 0-50±0.05 

-^ T = 1.0±0.1 ( — ) , (2) 

km s 1 \pc J 

gives S ~ 225 pc (C O line widths are 15kms _1 ; 
IWiklind & Combe^ lTggS). suggesting a cloud much larger 
than the 650 m as 2 solid angle subte nded by image A at 
2.9 GHz. While ISolomon et al.l ( fl987l) measure the relation 
for clouds with S < 50 pc and line widths < 8 km s" 1 , this 
is a strong indication that the cloud is q uite large, increasin g 
the probability of a high covering factor. iMulleret all ( 120071) . 
with high spectral resolution HCO + (2-1) observations, re- 



solve the line into four ~4.5 km s components, suggesting 
a clumpy medium of moderately sized clouds arrayed in front 
of image A, allowing the covering factor to remain constant 
or even increase with decreasing observation frequency (and 
increasing il A ). While the HCO + observations reveal that it 
is not a single large absorbing cloud, the presence of multiple 
clouds still allows for high covering factors. 

Agreement (lcr) between the observed optical depths and 
our LVG model (see Section l4~3l requires f^g > 0.2 for /14 5 = 
0.6 and /4 s > 0.3 for /14.5 = 1.0. Matching optical depths 
requires 0.5 < < 1 .0 for / 145 = 0.6 and / 4 8 = 0.5 for / 14 5 = 
1 .0. Although likely, / 4 g < /14 5 cannot be set a priori due to 
the clumpy absorbing medium. In subsequent analysis, we 
discuss results spanning the range of viable covering factors. 

4.3. Radiative Transfer Modeling 

Formaldehyde is susceptible to non-LTE excitation in either 
direction (maser inversion or anti-inversion), so it is generally 
inappropriate to assume a single line excitation temperature 
for the molecule. It is generally also not valid to attribute 
a single excitation temperature to all of the centimeter AJ = 
lines, since all line excitation temperatures can differ. We 
employ the still valid assumption of statistical equilibrium in 
these H2CO lines but allow all line excitation temperatures to 
float. 

We used an LVG radiative transfer model to find rela- 
tive level populations up to the 40th rotational energy state 
(IO37) of o-H 2 CO and the 41st (7 43 ) of p-H 2 CO. Assum- 
ing statistical equilibrium, the model computes optical depths 
and excitation temperatures for a given formaldehyde den- 
sity 77(H>CO) and a hydrogen density «(H2). The model as- 
sumes an isothermal and constant density spherically sym- 
metric cloud with large-scale turbulence or gravitational col- 
lapse producing large velo city gradients. Follow ing the ob- 
served velocity gradient of ISolomon et al.l (119871) . we set the 
velocity gradient to 1.0 km s" 1 pc " 1 . A^(o-H?CO) val ues scale 
inversely with the gradient dv/dr dWang et al .1120041) : 

W(o-H 2 CO) = [n(o-H 2 CO)/ (dv/dr)] x 3.08 x 10 18 x Av 

= [n(H 2 )-X /(dv/dr)] x 3.08 x 10 18 x Av, (3) 

where Av is the line width (FWHM), «(o-H2CO) is the num- 
ber density, and X is the abundance of 0-H2CO (not total 
H2CO) relative to H2. Since they are found via ratios of 
optical depths, which are proportional to A^(o-H2CO), n(H2) 
and T ex results are independent of dv/d r. For a more thor- 
ough discussion of t he LV G method, see lSoboievI d 19601) and 
iGoldreich & Kwanl ( fl974T) . 

Our model includes non-LTE collisional excitation of 
H2CO b y H2. We use the He collisional cross sections calcu- 
lated bv lGreenldl991l) as a proxy for the aspherical H2; Green 
notes that the substitution of He for H2 could cause errors up 
to 50% in individual collisional excitation rates between lev- 
els, but that this is likely to translate to overall errors in level 
populations of order 20%. Using pressure-broadened H2- 
H2CO and He-H2 CO laboratory measuremen ts of the mil- 
limeter transitions, iMengel & De Lucial d2000l) find the He- 
H2CO collision rates to be "very good," while the collisional 
cross sections of H2-H2CO collisions are erroneous by up to 
a factor of 2 from the scaled collision rates (yj Hh 2 I ^He) with 
no apparent pattern to the errors. 

We evaluated the effect on excitation temperatures of such 
an error in collision rates by randomly altering collision rates 
with Gaussian variation with a mean of zero and a standard 
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Figure 3. Contours trace excitation temperatures of the 1 iq— 1 1 1 (left) and 2 1 [ —2 1 2 (right) transitions as a function of «(H2) and jV(o-H2CO)/(km s pc~' )~' as 
calculated with our LVG model. Overplotted white lines represent the observed peak optical depth of the 1 10— 1 1 1 line (^horizontal) and the ratio of observed 
peak optical depths of the 1 10— In and 7.\^—1\i lines (^vertical) with 68% (1<t) confidence intervals. 7cmb = 4.6 K at z = 0.68, which we set as the background 
radiation temperatu re, and T^j, = 55 K I Henkel et al. 2005j). The black line is the 4.6 K contour. CS observations and modeling firmly limit niHj) < 2 X 10 4 cur 3 
IHenkel et al. 2005) — the range of these plots exceed that limit to show the behavior of H2CO as a function of «(H2) — and 7" cx < Tqmb f° r all points on the 
map below that limit. Column density (y-axis) is in units of cm -2 ( km s" 1 pc" 1 ) , and number density (.r-axis) is in units of cnT 3 . It is assumed throughout this 
paper that the velocity gradient dv/dr is 1 km s pc~' . The white at high «(H2) represents a ceiling to the contours at 12 K and not a physical plateau. The upper 
right corner of the 1 10— 1 1 1 plot reaches ~55 K, the kinetic temperature of the gas. Overplotted white lines on the 1 10— 1 1 1 plot represent the solution if covering 
factors are /4 s = 0.5 and /14 5 = 0.8; the central region is excluded, while a cloud with negligible collisions — \og(n[H2\) = cm~ 3 — but unrealistic X is not 
excluded by the ratio of optical depths at the 68% confidence level if priors on X are neglected. The similar lines on the 2\\—1\i plot represent the constraints 
placed by observations if covering factors are f$ g =0.3 and /14 5 = 1 .0; the central region is within the 68% confidence region although there is no solution in the 
model that exactly matches the observed lines. Differences in solutions shown on the two plots are due entirely to the selection of covering factors. 



deviation of 5 0%. T h is ove restimates the "up to 50%" er- 
ror quoted_by_(Gree^J 199l|) but is close to that measured 
by iMengel & De Lucia! (12000). For each of 25 sets of "er- 
roneous" collisional rates, we generated the parameter space 
shown in Figure [3] The effect of the errors on line excita- 
tion temperatures is s mall at low densities, where collisions 
are infrequent (Section l4~4l ); at moderate densities, the rms of 
the excitation temperatures of data sets is comparable to the 
mean excitation temperatures. For the density measurements 
reported herein (see below), the mean excitation temperatures 
of our erroneous data sets are 0.85 K for the 1 10 — 1 1 1 transi- 
tion and 1.4—1.7 K for the 2\\—2\2 transition; rms values are 
1.0 and 1.0-1.5, respectively. Taking an extreme scena rio and 
unifor mly doubling the collisional excitation rates of iGreenl 
( fl99lh increases the effect of collisions at a given density, ap- 
proximately halving the «(H2) values (0.3 dex) on the x-axis 
in Figure [3] 

The model used to determine excitation temperatures was 
generated with an ambient radiation temperature of Tcmb = 
2.73(1 +z) = 4.60 K and a kinetic temperature of 55 K. With 
a parameter space spanning 1 cm" 3 < «(H2) < 10 7 cm" 3 
and 3 x 10 12 cm" 2 (km s" 1 pc" 1 )" 1 < N(o-H 2 CO)/ (dv/dr) < 
3 x 10 16 cm" 2 (km s" 1 pc" 1 )" 1 with a sampling interval of 0.05 
dex, we compared the observed peak optical depths to those 
generated with our LVG model. The model is weakly de- 
pendent on kinetic temperature: changes in modeled optical 
depth are less than observational error for T% = 55 ± 20 K, and 
we present resu lts only for = 55 K, as measured via NH3 
observations bv lHenkel et al.l (l2005h. 

Of the four spectral components iMuller et al.l (|2007) iden- 
tify in their high resolution HCO + spectrum, the strongest ab- 
sorption occurs at the same velocity as both the H2CO 2n— 2u 
line and the strongest component of the 1 10 — 1 1 1 profile. The 
ratio of the strongest 1 10 — 1 1 1 absorption component and the 
2 1 1 —2 1 2 line allows in situ densitometric measurements: the 
optical depth of one transition gives A^(o-H2CO), and the ra- 
tio of two ^-doublets gives n(H2). Thus, with the two lines 
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Figure 4. Ratio of the optical depth of the 1 io~l 1 1 transition to that of the 
2 1 1 — 2 1 2 transition as a function of n(H2). The modeled optical depths assume 
covering factors of unity and a constant H2CO abundance relative to H2 of 
X = 10~ 9 . The two-line densitometry method cannot produce unique results 
in regions where a single value of the optical depth ratio corresponds to two 
values of ;i(H2) unless the density is constrained by priors or ancillary infor- 
mation. The non-unique region is enlarged if X is not assumed to be constant, 
as can be seen in the two viable density regions plotted in the left image of 
FigurefS] where the two possible densities are n(H2) ~ 6 X 10 4 cm' 3 (with 
X S3 10 -1 1 ) and n(H 2 ) ss 60 cnT 3 (with XRi3x 10~ 8 ). If priors can discrim- 
inate between ft(Hi) < 2.5 X 10 3 cnT 3 and ^(tb) > 2.5 X 10 3 cm -3 , solutions 
are unique and tbCO is a powerful densitometer. 



one can define a region of the modeled parameter space by 
matching observed optical depths to modeled ones. When 
n(H2) <2x 10 4 , the ratio of lines can result in bifurcated solu- 
tions for n(H2) when using the two-line densitometry method 
employed herein — one at moderate n(H2) and one at low 
n(H2), as is shown in Figure|4] 

To break the degeneracy of solutions using observations 
of other molecules and abundance arguments, we divide the 



Formaldehyde Anti-Inversion in B02 18+357 



7 



17 



\ 
> 

o 

w 

I 

I 

o 14 



60 
O 



II 



III 



/■ / A / / 



' ! / / ! 

/ / ' t/ ' 

/ / i / /i / 

/ / i / / i ' 

/ / i . X / i ' / 

l/ * *f / / 

' / \ / / \ / y 

/ A / * / ' 



/ 



13 E- ' 

/ 



12 3 4 5 6 

log(n[Hj) 

Figure 5. Diagonal dashed lines mark the abundance of 0-H2CO relative to 
H2 (see Equation|5J; ve rtical dot-dashed lines separate the regions, labeled at 
top, discussed in Section l4~3l and Table[3] Typical Galactic abundances are of 
order 1(T 9 , while measurements of Af(H 2 CO ) IJethava et al.120071) and N(H 2 ) 
(Combes & Wiklind 1995; Gerin et al 1997) place the total H 2 CO abundance 
between 2.5 X 10^ and 10~ 10 in B02 18+357; the ortho/para ratio, which 
Jethava et al. (2007) measure as ~2.8, drops the abundance of 0-H2CO by 
~25% relative to the total abundance of EE>CO. 



possible solutions into three regions separated at «(H.2) = 
2 x 10 2 cm" 3 and n(H2) = 2 x 10 4 cm" 3 , as shown in Fig- 
ure [5] and described in Table [3] The Region I/1I boundary at 
n (H2) = 2 x 10 2 is an arbitrary division between low densities, 
where T ex rj Tcmb f° r both the I10-I11 and 2n—2\2 transi- 
tions, and moderate densities, where T ex < Tcmb- The Re- 
gion Will boundary at nQ&i) = 2 x 10 4 cm" 3 is t he firm upper 
limit placed by CS observations and modeling dHenkel et alj 
l2005h . Ranges of abundances in each region are calculated 
for all covering factors for which equilibrium model solu- 
tions exist, determining A^o-rTCO) and n{Y\2) for all cover- 
ing factors within these ranges, and converting column den- 
sities to number densities of 0-H2CO with Equation [3j max- 
imum and minimum abundance values (X, Equation |3]l listed 
in Table |3] are the maximum and minimum values of n(o- 
H2CO)/«(H2) in each region. Due to the high abundance 
(log[X] > -8.2), the low-density Region I is unlikely: typical 
Galactic abundances have a small dispersion around a few x 
10" 9 (e.g., lEvans et al.1 H975t iDickel. Goss. & Rotl [T987h . 



Measurements of jV(H 2 CO) (iJethava et alj|2007h and N(H 2 ) 
dCombes & WikUndl fl995l : iGerin et all 11997b in B02 18+357 
place the total H2CO abundance between 2.5 x 10~ 9 and 
10~ 1Q in B02 18+357; the ortho/para ratio, which lJethava et alj 
d2007h measure as 2.8 in their best model of the lines they 
observed, drops the abundance of 0-H2CO by ^25% from the 
total H2CO abundance. These constraints agree with the value 
of «fH 2 ) w 3.2 x 10 3 cm" 3 estimated from NH 3 dHenkel et alj 
2005, no confidence interval given). 

Fitting observed optical depths to our LVG model within 
the above constraints and forcing /14.5 > 0.6 gives 2 x 
10 3 cm" 3 < n(H2) < 1 x 10 4 cm" 3 (lcr) for the covering fac- 
tors discussed above. Agreement between observations and 
modeling requires /4.8//14.5 > 1/3, while the ratio of solid 
angles in image A is ^4.8/^14.5 ~ 20. Since the covering 
factors are poorly constrained, results for n(H2) and T ex can- 
not have proper confidence intervals. However, our reported 
ranges include all reasonable covering factors (/14.5 > 0.6, 
f4.&/ /14.5 1/3) and are thus robust. We list n(H.2) and T ex 



Table 3 

Formaldehyde Radiative Transfer Model Regions: H2 Number Density and 
H2CO Abundance 



Region 


Min. log(»(H 2 )) 


Max. log(n(H 2 )) 


Min. log(X) 


Max. log(Z) 


I 


0.0 


2.3 


-8.2 


-5.6 


II 


2.3 


4.3 


-10.2 


-8.2 


III 


4.3 


5.5 


-10.3 


-10.2 



Note. — Molecular hydrogen densities ^(Eb), with units of cm , and o-H 2 CO 
abunda nces X (Equation |5J at the boundaries between the regions discussed in 
Section l431 and shown in Figure|5] The Region I/II boundary is an arbitrary di- 
vision between low densities, where T ex ~ Tcmb . ™d moderate densities, where 
T ex ^ Tcmb ■ The R egion II/III boundary is the upper limit placed by CS obser- 
vations and modeling Henkel et al. (2005). Ranges of abundances are calculated 
by determining A r (o-H 2 CO) and ;i(H 2 ) for all covering factors with equilibrium 
model solutions within these ranges, and converting column densities to number 
densities of 0-H2CO with EquationfJ] maximum and minimum X values are the 
maximum and minimum values of ;i(H2CO)A?(H 2 ) in each region. The typical 
H2CO abundance relative to Eb in Galactic molecular clouds is ~10~ 9 , and the 
ortho/para ratio ranges from 1 to 3. In B02 18+357, the total H2CO abundance is 
1 —25 X 10~ 10 , assuming an ortho/para ratio of 2.8. 




3 4 
log(n[H 2 ]) 

Figure 6. Excitation temperatures of the lowest five ^-doublet transitions as 
a function of h(H 2 ) for W(o-H 2 CO) = 10 14 cm 4 , as calculated with our LVG 
model. This figure is a horizontal slice through Figuref3]with additional K- 
doublet lines added. The divergence of the ^-doublet excitation temperatures 
as a function of 71(1+2) makes the ratio of the any pair of lines a sensitive den- 
sitometer. However, since the ratios converge to unity at low >i(H 2 ) (as T cx 
approaches Tcmb) ™d cross it again at fl(H2) ~ 10 s (as all excitation temper- 
atures go to Tki n at high densities), the ratios cannot provide a unique global 
solution unless the covering factor is known a priori or other information, 
such as obse rvations of other molecules, is available to break the degeneracy 
(Sectionl4~3l. 



for various covering factors in Tables |J-|6| T e 
physical pairs of covering factors. 



< T C mb for all 



4.4. Column Density and Excitation Temperatures 

We find 2.5 x 10 13 cm" 2 < jV(o-H 2 CO) < 8.9 x 10 13 cm" 2 
(assuming a velocity gradient of 1 km s" 1 pc" 1 ), with the 
range determined by the viable covering factors (/14.5 > 0.6, 
/4.8//14.5 ^1/3) and constraints on molecular hydrogen den- 
sities (2 x 10 3 cm" 3 < «(H2) < 1 x 10 4 cm" 3 ). Since these re- 
sults cover all reasonable covering factors, the limits are hard 
and conservative. The total H2CO column dens ity, assuming 
0-H2C O/P-H2CO = 2.8 from the best models of IJethava et alj 
( 120071) , which is determined primarily from millimeter obser- 
vations where covering factors are unity and not a concern, 
is 3.4 x 10 13 — 1.2 x 10 14 cm" 2 . The column density scales 
inversely with the assumed velocity gradient. 
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Table 4 

Model-Derived Number and Column Densities For 
Various Covering Factors 



/4.8 


fl4.5 


log(n(H 2 )) 
(cm" 3 ) 


log(JV(o 
(cm -2 


-H 2 CO)/(dv/dr)) 
(kms" 1 pc-'r 1 ) 


0.50 

0.50 

0.30 
1.00 


0.60 

0.80 

1.00 
1.00 


5 05 rt 15 

0- 03^ 
4.80^j° 

1 80+ 045 

1 - -1.80 

3 53+0- 37 
5 15 mi5 




1 J - -0.10 
i 4 iA+0.00 

14 - lu -o.io 
n 4 5 +o.io 

1J - -0.15 
1 3 60+ 005 

13 30 +0 1 ° 



Note. — Confidence intervals are set at \a (68%) 
but are non-Gaussian. Where more than one solution ex- 
ists within lcr observational error, all solutions are given. 
An eiTor of '0.00' indicates that the modeled density 
did not change significantly across the region of density 
space defined by fitting observed optical depths to those 
derived with our LVG model. Number densities above 
2 X 10 4 cm -3 contradict the firm upper limit obtained from 
LVG models of CS observations (Henkel et al. 20051). 

Regardless of covering factors (although tables herein only 
list several, the entire range of reasonable covering factors 
was evaluated), Figure [3] shows that the excitation tempera- 
tures of both observed lines are below Tcmb within the con- 
straints placed by CS, which require n(H2) <2x 10 4 cm" 3 . 
Figures [3] and [6] show clearly that, at z = 0.68, the H2CO cen- 
timeter lines will be anti-inverted. The non-LTE excitation is 
density dependent: 

1. n(H2) < 10 2 cm" 3 : H2CO excitation is dominated by 
CMB photons. All excitation temperatures are held 
near to the microwave background and the system is 
nearly in radiative equilibrium with the CMB because 
collisions with H2 are too infrequent to pump the H2CO 
^-doublet population into anti-inversion. 

2. 10 2 cm" 3 < n(H 2 ) < 3 x 10 5 cm" 3 : Collisions with H 2 
anti-invert the H2CO ^-doublet population. The mil- 
limeter transitions, however, have T ex > Tcmb- 

3. «(H2) >3x 10 5 cm" 3 : High densities increase collision 
rates and drive H2CO toward thermal equilibrium with 
H2, with a kinetic temperature of 55 K. 

Figure |6l a horizontal slice through Figure [3] shows a 'IT 
shape of excitation temperatures. The two possible values 
of T ex and n{Hi) for a given ratio of obser ved optical depths 
cause the degeneracy discussed in Section 14.31 and shown in 
Figure|4] The n(H2)-dependence of centimeter excitation tem- 
peratures makes H2CO a sensitive in situ densitometer if the 
degeneracy can be broken with a priori knowledge of physical 
conditions such as H2CO abundance or limits on «(H2). 

5. DISCUSSION 

5.1. Comparison to Excitation in Galactic Molecular Clouds 

Galactic dark clouds show a decrement in the excitation 
temperature of the lip— I n line relative to th e CMB of T ex - 
Tcmb - -0.7 to -0.4 K (lEvans et al.1 [l975h . An LVG ra- 
diative transport model applied to Galactic giant molecular 
cloud s derives comparab le decrements, ranging from -1 to 
K dHenkel et al.1 1 19801) . The cloud studied in this paper 
shows a much larger decrement, with T ex -7cMB ^5 _ 3.6 K (for 



2 x 10 3 cm" 3 < «(H 2 ) < 1 x 10 4 cm" 3 ). For the 2 n -2i 2 transi- 
tion, the decrement is ~ -2.6 to -3. 1 K. The magnified decre- 
ment is due to the coll isional dependence of excitation tem- 
peratures (Section 14.41 1: CMB excitation is only significant 
in diffuse clouds with «(H2) < 10 2 cm" 3 , so for higher density 
clouds the excitation temperature is nearly constant regardless 
of the CMB temperature. Increasing the background temper- 
ature thus magnifies the T ex - Tcmb decrement proportionally 
to(l+z). 

Because T ex - Tcmb < for both transitions, the cloud 
would be observable in absorption against the CMB in both 
lines without the background quasar - if the cloud were in a 
different place in the lensing galaxy rather than aligned with 
image A or in a galaxy with no background continuum source, 
for instance. In the case where the cloud is matched to the 
telescope beam (no beam dilution) at z = 0.68 and the cov- 
ering factor is unity, the observed 1 10 — 1 1 1 line temperature 
would be 

7ii„e = (l-e- T )(T ex -T CM B)/(l+z) » "36 mK. (4) 

For a non-unity covering factor and small optical depth, as 
seen in B02 18+357, 7ii ne is approximately inversely propor- 
tional to the covering factor. Such a cloud illuminated by the 
CMB will be observable with future instruments such as the 
EVLA. Since T ex is small and nearly constant with redshift 
and Tcmb scales as 1 +z, this raises the possibility of survey- 
ing molecular gas in galaxies regardless of distance or chance 
alignment with sources of background illumination. 

Models indicate that the anti-inversion persists in the K a = 1 
/^-doublets at least up to the 615-616 0-H2CO transition for all 
reasonable covering factors (which determine the measured 
molecular hydrogen density), and for some covering factors 
the anti-inversion persists into the transition (Table 

|5]). Although the highest energy lines may be difficult to de- 
tect, all of the anti-inverted transitions should in principle be 
observable in absorption against the CMB. 

^-doublet transitions with K a = 2 (p-ITCO) and K a = 3 
(0-H2CO) are even more strongly anti-inverted than those in 
the K a = 1 ladder, but relative to the K a = 1 /if-doublets the 
higher energy of the K a > 1 states (the lowest K a = 2 K- 
doublet level, 221 at 57.6 K, has approximately the same en- 
ergy as the 5 15 level at 62.5 K) generally makes them more 
optically thin and difficult to observe given the low temper- 
atures of molecular clouds. Our models indicate that the 
220-221 transition (71.1 MHz rest frequency), the lowest K- 
doublet of the K a = 2 ladder, has an excitation temperature 
less than 10 mK in B02 18+357. Due to the relatively warm 
temperature of the gas observed in B02 18+357 (55 K), the 
optical depth of the transition is predicted to be lower than 
that of the observed 1 10 — 1 1 1 0-H2CO line by only a factor of 
r^5 and higher by a factor of ^20 than would be expected 
in an otherwise equivalent cloud with a kineti c temperature 
of 15 K (assuming an ortho/para ratio of 2.8; iJethava et alj 
120071) . The next lowest ^-doublet of the K a = 2 ladder, the 
32i _ 322 transition (355.6 MHz rest), is predicted to have an 
optical depth that is an order of magnitude lower than that 
of the 220-221 line. All other K a > 1 /^-doublets are weaker 
by yet another order of magnitude or more than the 3ti— 322 
line. In the hypothetical case of a similar gas-rich galaxy il- 
luminated only by the CMB at z = 0.68 (the flux from the 
quasar illuminating B02 18+357 is not well known at meter 
wavelengths), the observed line temperatures for the lowest 
K a = 2 A'-doublets are predicted to be between -4 and -11 
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Table 5 

Modeled A^-doublet ( AJ = 0) 0-H2CO Line Excitation Temperatures 



As /14.5 log(n(H 2 )) 1 10-1 11 
(cnT 3 ) 4.8 GHz 



2ii~2i2 3i2— 3i3 4 13 -4 14 5 14 — 5 15 6 15 — 6 16 7 16 -7 17 
14.5 GHz 29.0 GHz 48.3 GHz 72.4 GHz 101.3 GHz 135.0 GHz 



0.50 0.60 5.05+ 



u - u:) -0.03 



0.50 0.80 



0.30 
1.00 



1.00 
1.00 



4.80!°; 2 ° 

3 53+°- 37 
-0.38 

5 15+ 015 



97+3-0 
u - ''-0.09 
4 57+O.OO 
^■•"-1.23 

0-84^ 
3 24 + 1 - 33 



1 57+z.yu 
-0.03 
4 5:0+0.00 

+2.95 
-0.00 



1.53 



4 qt +0.56 

1 75+2.68 
i./J_0 10 

1 60+ 291 
i.ou_ 003 



1 98 +/ 33 
4 50+0.00 

^■ J? -0.28 

1 98+ 2 - 57 

4 30+ - 29 

2 51+ 2 - 01 

-0.20 

20 1+2.55 
z - ul -0.03 



2 34+ 2 - 22 

J -0.04 
4 60+ 000 

9 90+2.27 

^■ i7 -o.oo 

^' -0.30 

9 »q+1.66 
-0.27 
9 aq+2.18 



9 Q4+1- 61 

4 fiO+o 00 

^■°"-0.20 

9 70+1. 77 
z -''-0.07 
4 ao+0.21 
^•■"-0.31 

2 - 87 ^:?6 

3.03^ 



3 63+"- 10 
J - -0.14 

4 S 8+0-00 

14^+1.00 
-0.12 
"3 09+O.66 

3 06+ 129 
^■^-o.oo 

3 76+ - 71 



4 77+0.20 
^■"-0.79 
441+O.OO 
^•^ 1 -0.48 
4 ci+0.20 
^■ ; '-'-0.53 
t Q3+0.48 
J ' 7J -0.03 

T. Q3+0.07 

J? -0.04 

4 92+0- 25 
^■^''-0.90 



Note. — Excitation temperatures for selected pairs of covering factors fa g and /14.5; boldfaced temperatures have 
T ex > Tcmb (Tcmb = 4.60 K) within the listed 68% (non-Gaussian) confidence intervals. Frequencies are rest fre- 
quencies, and all excitation temperatures have units of K. Regions listed with errors of '0.00' do not change excitation 
temperature significantly across the region of density space defined by fitting observed optical depths to those derived 
with our LVG model. Multiple excitation temperatures for a given pair of covering factors correspond to solutions in 
Table|4]where a pair of optical depths did not produce a unique solution. Low density solutions MH2) < 10 2 cm* 5 ) have 
T cx = Tcmb (z = 0.68466) = 4.60 K because collisions are insufficient to be a significant excitation mechanism. Of the 
^-doublet transitions, only the lio - In and 2n— 2(2 lines have been observed in B0218+357 (see Tables- 



Table 6 

Modeled AJ = 1 H2CO Line Excitation Temperatures 



As /14.5 log(n(H 2 )) 2 12 -1„(0) 2„-l 10 (O) l m -0oo (P) 2 02 -loi (P) 
(cm 4 ) 140.8 GHz 150.5 GHz 72.8 GHz 145.6 GHz 



0.50 


0.60 


5 05+0- 15 
J - UJ -0.15 


7 7^+1-45 
'•' -3.13 


6 70+0- 89 
-2.10 


,7 40+14.90 
11 -^"-12.80 


6 74+0-84 
D -' 4 -2.14 






0.Q£j;g 


4 60+° 00 
^■ ou -o.oo 


4 60+0-0° 
^■ ou -o.oo 


^■ DU -o.oo 


4 60+0- 00 
^ D -0.00 


0.50 


0.80 


4.80+o : 2 o 


6 31+1 06 


5 7O+0.69 
5./8_j 18 


q 77+5.23 


5.84+ ; 67 








4 60+° 00 
^■ ou -o.oo 


4.60!°;«° 


4-«£gffi 


4 60+0- 00 
^ DU -0.00 


0.30 


1.00 


3 53+0- 37 


4-73^;?f 




4.88+0;'! 


4 70+0- 05 
^■ /u -o.io 


1.00 


1.00 


5 15+0- 15 
J - 1J -o.io 


8.63^ 




25-20^;™ 


7 97+I.O8 
'• z '-2.67 



Note. — Similar to Table [3] but for AJ = 1 transitions of both o-H 2 CO (O) and 
p-H 2 CO (P). All excitation temperature s have units of K. All transitions in this ta- 
ble have been observed in B02 18+357 I Jethava et al. 2007). Low density solutions 
(n(H 2 ) < 10 2 cm" 3 ) have T ex = Tcmb (z — 0.68466) — 4.60 K because collisions are 
insufficient to be a significant excitation mechanism. 



mK for the 2 20-221 line and -0.3 and -1 mK for the 3ti— 322 
line (assuming YV(H2CO) = 5 x 10 13 cm" 2 and an ortho/para 
ratio of 2.8, with the range set by our measured range of val- 
ues for the molecular hydrogen density in the cloud). While 
these low-frequency lines will be dwarfed by the foreground 
Galactic continuum and difficult to observe for their own sake, 
this modeling indicates that H2CO K a = 2 absorption lines in 
gas-rich galaxies present potential contaminants to Epoch of 
Reionization studies. 

5.2. Comparison to Previous Results 

Uethava et all (120071) measure M0-H9CO) = 3.8 x 10° cm" 2 
and «(H 2 ) < 10 3 cm" 3 (with a best fit of n(H 2 ) = 200 cm" 3 ) 
with their n on-LTE LVG analysis of six H2CO lines in 
B02 18+357. iMenten & Reidl (119961) . assuming LTE, mea- 
sure YV(o-H2CO) = 1.2 x 10 13 cm" 2 from their observation 
the 2\\-2\2 line, although they estimate that anti-inversion 
of the ^-doublet would reduce their res ult by 50%, well be- 
low the column densities reported by iJethava et al.l (120071) 
and our range of 2.5 x 10 13 -8.9 x 10 13 cm" 2 (for dv/dr = 
1 km s" 1 pc" 1 ). 

Our 1 io~l 11 spectrum has a significantly higher signal to 
noise ratio and better spectral resolution than previous obser- 



vations, and reliance solely on it and our new 2n-2!2 data 
to determine densities precludes high-frequency atmospheric 
variability while magnifying the effect of uncertain cover- 
ing factors. Since the centimeter ^-doublet lines measure 
«(H2) with greater p recision than the mil limeter lines can, the 
trade-off makes the ' Jethava et al.l (120071) result less accurate 
for determining n(H-2) and T ex (which they do not address) 
but more reliable for measuring AMo-HbCO), since they rely 
heavily on nearly thermal millimeter transitions which have 
covering factors of unity. The optimal combination of the 
two methods uses the centimeter lines to constrain n(H-2) (and 
T ex , since it is closely related) and the millimeter lines to con- 
strain A^o-LLCO). T ex and n(H.2) are nearly independent of 
YV(o-H2CO) at the column densities observed in this source 
(Figure |3J, 

Our grid-search method is a more comprehensive analy- 
sis of the »(H?)-/y( o-H?CO) density space than that used by 
IJethava et al.l (120071) . who compared their observations only to 
LVG models with «(H 2 ) = 0, 200, 500, 1000, and 3000 cm" 3 
to determin e a best-fit (y 2 ) valu e. The small values of «(U-2) 
sampled by iJethavae t al. (20 07) a lso neglect the high-n(H.2) 
solutions discussed in Section POI essentially assuming their 
conclusion that the absorbing cloud is diffuse and low-density 
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and not considering the alternative higher density regime. 

6. CONCLUSIONS 

With new measurements of the lio— In an d 2n— 2i2 transi- 
tions of H2CO in B02 18+357, we have measured the molecu- 
lar hydrogen number density and H2CO column density. Us- 
ing an LVG model of the physical conditions of the absorbing 
cloud in the gravitational lens and assuming a velocity gra- 
dient of 1 km s" 1 pc" 1 , we measure /V(o-H2CO) to be in the 
range 2.5 x 10 13 -8.9 x 10 13 cm" 2 , where the uncertainty in 
values is caused by the uncertainty in covering factors (in the 
optically thin limit, the colum n density is i n versel y propor- 
tional to the covering factor). iJethava et al.l (120071) measure 
an ortho/para ratio of 2.0-3.0 (with a best value of 2.8) and 
a second component to the profile that adds ~15% to the to- 
tal column density; these two contributions increase the to- 
tal column Af(H2CO) by ~55% above our measurement of 
TVYo-FLCO) in the primary absorption component. Since col- 
lisions with hydrogen are the dominant excitation mechanism, 
we are also able to measure n(Hi) = 2 x 10 3 - 1 x 10 4 cm" 3 , 
which again spans the range of viable covering factors. 

Modeled excitation temperatures show both observed cen- 
timeter lines to be anti-inverted relative to the 4.6 K CMB at 
z = 0.68, with T„(lio-ln) w 1 K and T ex (2 n -2i 2 ) < 2 K. 
Excitation temperatures of the millimeter lines are >4.60 K, 
as expected. These results show that the dominant absorbing 
cloud in the B02 18+357 lens is analogous to Galactic molec- 
ular clouds that show anti-inversion of the /^-doublet transi- 
tions, and the centimeter H2CO lines could be observed in 
absorption against the CMB if the background source was re- 
moved or if the cloud was viewed along a different line of 
sight. 

Since T ex is determined by collisions in dense clouds, the 
source-frame decrement T ex - Tcmb is magnified when Tcmb 
is increased (n.b., at highredshift). The anti-inverted centime- 
ter transitions thus could provide absorption lines against the 
CMB that would trace molecular gas independently of dis- 
tance when observed with a telescope capable of resolving 
the star-forming regions in galaxies. Future work targeting 
additional objects will quantify the change in T cx -Tcmb with 
redshift and test these results. 

We are indebted to C. Henkel for providing the LVG code 
used in this analysis and to the anonymous referee who sug- 
gested that we model the excitation temperatures of K a > 1 
/^-doublets and checked our calculations with an independent 
radiative transfer code. 
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